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ABSTRACT 

Aimed at ionization balance and spectral analysis of UV and X-ray sources, 
we present self-consistent sets of photoionization cross sections, recombination 
cross sections, and rate coefficients for Li- like O VI and He- like O VII. 
Relativistic fine structure is considered through the Breit-Pauli R-matrix 
(BPRM) method in the close coupling approximation, implementing the unified 
treatment for total electron-ion recombination subsuming both radiative and 
di-electronic recombination processes. Self-consistency is ensured by using an 
identical wavefunction expansion for the inverse processes of photoionization 
and photo-recombination. Radiation damping of resonances, important for 
H-like and He-like core ions, is included. Compared to previous LS coupling 
results without radiative decay of low-n n < 10) resonances, the presents results 
show significant reduction in O VI recombination rates at high temperatures. 
In addition to the total rates, level-specific photoionization cross sections and 
recombination rates are presented for all fine structure levels n {£ SLJ) up to n 
< 10, to enable accurate computation of recombination-cascade matrices and 
spectral formation of prominent UV and X-ray lines such as the 1032,1038 A 
doublet of O VI, and the 'triplet' forbidden, inter combination, and resonance 
X-ray lines of O VII at 22.1, 21.8, and 21.6 A respectively. Altogether, atomic 
parameters for 98 levels of O VI and 116 fine structure levels of O VII are 
theoretically computed. These data should provide a reasonably complete set of 
photoionization and recombination rates in coUisional or radiative equilibrium. 

Subject headings: atomic data — atomic processes — photoionization, 
dielectronic recombination, unified electron-ion recombination — X-rays: 
general — line:formation 
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1. INTRODUCTION 



Li-like and He-like Oxygen are among the most important atomic species in UV and 
X-ray plasma diagnostics of astrophysical sources. For example, there is considerable 
current interest in O VI and O VII as a possible reservoir of 'missing baryons' in the 
warm/hot intergalactic medium in the present universe (Cen et al. 2001, Fang and Canizares 
2000, Yoshikawa et al. 2003). In addition to global ionization balance models that require 
O VI/0 VII photoionization and recombination rates, specific spectroscopic diagnostics of 
both of these ions are related in an interesting but intricate manner. Outflows observed 
from the 'warm absorber' in active galactic nuclei exhibit absorption/emission lines in the 
X-ray from KLL resonant photoabsorption in O VI, and excitation/recombination of O VII 
lines. This phenomenon was investigated theoretically by Pradhan (2000) who pointed out 
that the observed spectra of the well known AGN NGC 5548 shows these features (Kaastra 
et al. 2000). Recently Arav et al. (2003) have further pursued this suggestion to infer large 
discrepancies in the UV and X-ray O VI column densities in the outflow from NGC 5548. 
Ionization balance and spectral formation of O VI and O VII are crucial factors in all such 
studies. 

The inverse atomic processes of photoionization and clcctron-ion photorecombination 
of Li-like and He-like ions are also of particular interest in X-ray astronomy to 
analyze new observations by space based observatories such as the Chandra X-ray 
Observatory and and XMM-Newton, at photon energies and temperatures prevalent 
in high-temperature sources such as AGB, supernova remnants, hot stellar coronae 
etc. (e.g. Canizares et al. 2000, in Proceedings of the NASA workshop on Atomic data 
needs in X-ray astronomy 2000). X-ray emission in the Ka complex comprising of the 
'triplet' feature of resonance (w), intercombination (x,y), and forbidden (z) lines in 
He-like ions, corresponding to the 4 transitions from the n — 2 levels to the ground 
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level Is^ (^So) < — ls2p(^P°),ls2p(^P°i),ls2s(^Si) respectively, yields valuable spectral 
diagnostics of temperature, density, ionization balance, and abundances in the plasma 
source. Calculations of the recombination- cascade contributions for these important lines 
requires accurate atomic parameters for fine structure levels up to fairly high n levels, as 
presented in this report. 

In the Present work we report new results for self-consistent data sets for photoionization 
and recombination for (e -|- O V) — > O VI and (e -|- O VI) — > O VII, using a unified 
theoretical method that takes account of both radiative and di-electronic recombination 
processes (RR and DR). Although treated separately in other previous treatments, both 
recombination processes, RR and DR, are inseparable in nature. A unified treatment 
is therefore not only theoretically and computationally more accurate, it is also more 
practically convenient for applications since a single recombination rate coefficent takes 
account of both RR and DR in an ab initio manner. For highly charged ions it is important 
to consider relativistic fine structure explicitly in the theoretical formulation, in addition to 
the electron correlation effects. Previous results for O VI (Nahar 1998, 1999) were obtained 
in a non-relativistic LS coupling approximation; relativistic effects were partially included 
for O VII. To further improve the accuracy, as well as to provide extensive sets of data for 
astrophysical models, we have re-calculated all photoionization cross sections, electron-ion 
recombination cross sections, and rate coefficients using the relativistic Breit-Pauli R-matrix 
(BPRM) method. Similar calculations have previously been reported for He- and Li-hke 
ions C V and C IV (Nahar et al. 2000), and Fe XXV and Fe XXV (Nahar et al. 2001). 

2. THEORY 

The unified method of electron-ion recombination (Nahar and Pradhan 1992) yields 
self-consistent sets of atomic paramters for photoionization and recombination for atoms 
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and ions (e.g. Nahar and Pradhan 1997, Paper I of this series). Relativistic extension using 
the BPRM method is described in Zhang and Pradhan (1997), Pradhan and Zhang (1997) 
and Zhang et al. (1999). Photorecombination of an incident electron with the target ion 
may occur through (i) non- resonant, background continuum, or radiative recombination 



which is the inverse process of direct photoionization, or (ii) through a two-step 
recombination process via autoionizing resonances, i.e. dielectronic recombination (DR): 



where the quasi-bound doubly-excited autoionizing state may lead either to (i) 
autoionization, a radiation-less transition to a lower target state with the electron going 
into a continuum, or (ii) radiative stabilization to a recombined bound state via decay of 
the ion core (usually to the ground state) with the electron captured. 

The extension of the R-matrix method to electron-ion recombination calculations 
entails close coupling calculations for photoionization and electron-ion scattering. Identical 
eigenfunction expansion for the target (core) ion is employed for both processes, enabling 
inherently self-consistent results in an ab initio manner for a given ion. We consider 
photoionization from, and recombination into, the infinity of levels of the (e -|- ion) system. 
These are divided into two groups of bound levels: (A) with v < Vq and all possible fine 
structure Jtt symmetries, and (B) z/q < i/ < oo; where v is the effective quantum number 
relative to the target threshold (s). Photoionization and recombination calculations are 
carried out in detail for all group A levels, while group (B) levels are treated through 
quantum defect theory of DR within close the coupling approximation (Bell and Seaton 
1985, Nahar and Pradhan 1994). A generally valid approximation made in recombination 



(RR), 



e + X++ ^hu + X+, 
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to group (B) levels is that the background contribution is neghgible, and DR is is the 
dominant process in the region below the threshold of convergenece for high-n resonances. 

In the close coupling (CC) approximation the target ion (core) is represented by an 
A^-electron system. The total wavef unction, "^{E), of the (A^+1) electron- ion system of 
symmetry Jtt is represented in terms of an expansion of target eigenfunctions as: 

^iE)=AJ2x^0^ + J:C,^„ (3) 
i j 

where Xi is the target wavefunction in a specific level JjTTi and 9i is the wavefunction for 
the (A'"-|-l)-th electron in a channel labeled as SiLi{Ji)7:ik^ii{ Jtt)); kf being its incident 
kinetic energy. $j's are the correlation functions of the (A'"-|-l)-electron system that account 
for short range correlation and the orthogonality between the continuum and the bound 
orbitals. In relativistic BPRM calculations the set of SLn are recoupled for Jtt levels of (e 
+ ion)-system, followed by diagonalisation of the Hamiltonian, 

H^^.W = EW, (4) 

where the BP Hamiltonian is 

iirBP [7jV_K I TTinass , rrDar i rrso /r\ 

The first term, -ff^+i, is the nonrelativistic Hamiltonian, 

and the additional one-body terms are the mass correction term, H'^^^^^ = —^j^ipj, 
the Darwin term, iJ^ar _ ^^^^.^/'^(^l-y and the spin-orbit interaction term, 
— Zo? Yl,i ^li-Si., respectively. 

The positive and negative energy states (Eq. 4) define continuum or bound (e -|- ion) 
states, such that, E — k'^ > ^ iox continuum (scattering) channels and E — < 
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for bound states, where v is the effective quantum number relative to the core level. If 
£^ < then all continuum channels are 'closed' and the solutions represent bound states. 
The matrix element for the bound-free transition, < ^'sHDH^f >, can be obtained from 
the continuum wavefunction {^f) and the bound wavefunction {<^b)- D is the dipole 
operator, = Zli'^i, in length form and Dy = — 2^^^ Aj in velocity form, the sum being 
the number of electrons. 

Photoinization cross sections of all possible states with n < rimax ~ 10 - are obtained 
as in the CC approximation as in the Opacity Project (Scaton 1987, The Opacity Project 
1995, 1996), extended to include the relativistic effects under the Iron Project (Hummer et 
al. 1993). The photoionization cross section is obtained as 

a PI = -—ujS, (7) 

g 3c 

where g is the statistical weight factor of the bound state and S is the dipole line 
strength, S = | < ^^bHDH^j? > p. For highly charged H- and the He-like ions, the 
probability of radiation decay of an autoionizing state is usually comparable to that of 
autoionization as discussed in Nahar et al. (2000). With strong dipole allowed 2p — > Is 
and ls2p (^Pi) — > Is^ C^o) transitions autoionizing resonances are radiatively damped to 
a significant extent. The radiative damping effect of all near-threshold resonances, up to 
u < 10, is included using a resonance fitting procedure (Sakimoto et al. 1990, Pradhan and 
Zhang 1997, Zhang et al. 1999). 

The photo-recombination cross section, a^c, is related to photoionization cross section, 
CTpi, through principle of detailed balance (Milne relation) as 

<7Rc(ej = CTpi [6) 

in Rydberg units; a is the fine structure constant, e is the photoelectron energy, and / 
is the ionization potential. (Trc are computed from the photoionization cross sections 
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at a sufficiently large number of energies to delineate the non-resonant background and 
the autoionizing resonances, thereby representing both radiative and the dielectronic 
recombination (RR and DR) processes. In the unified treatment the photoionization 
cross sections, crpi, of a large number of \ow-n bound states - all possible states with 
n < rimax ~ 10 - are obtained as described above. It is assumed that the recombining ion is 
in the ground state, and recombination can take place into the ground or any of the excited 
recombined (e+ion) states. 

Recombination rate coefficients of individual levels are obtained by convolving 
recombination cross sections over Maxwellian electron distribution f{v) at a given 
temperature. 



Contribution of these \ow-n group (A) bound states to the total ctrc is obtained by summing 
over all low-n group (A) states, ctrc thus obtained from crpi, including the radiatively 
damped autoionizing resonances (Eq. (8)), corresponds to the total (DR+RR) unified 
recombination cross section. 

Recombination into the high-n states, nmax < n < oo, (Fig. 1 of Nahar & Pradhan 
1994) are included separately. To each excited threshold of the core, JiTTj, belongs an infinite 
series of (A^+l) -electron levels, JiTTiui, to which recombination can occur. For the high u 
levels DR dominates while the background RR is negligibly small. The contributions from 
these levels are added by calculating the DR collision strengths, fl^R, employing the precise 
theory of radiation damping by Bell and Seaton (1985; Nahar & Pradhan 1994). The 
recombination cross section, auc in Megabarns (Mb), is related to the collision strength, 
flRC, as 



where k'f is the incident electron energy in Rydbergs. As ctrc diverges at zero-photoelectron 




(9) 



<TRc{t ^ j){Mb) = 7^0«c(^, j)/(5i^-)(«o/l- x 10"''), 



(10) 
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energy, the total collision strength, Q, is used in the recombination rate calculations. 
Background photoionization cross sections for the high-n group (B) levels are computed 
hydrogenically, which is referred to as the "high-n top-up" (Nahar 1996). 

3. COMPUTATIONS 

Relativistic BPRM calculations in intermediate coupling are carried out in the close 
coupling approximation using the i?-matrix package of codes. These codes are extensions of 
the Opacity Project codes (Berrington et al. 1987) to include relativistic effects (Scott & 
Burke 1980, Scott & Taylor 1982, Berrington et al. 1995), as implemented under the Iron 
Project (Hummer et al. 1993). The calculations span several stages of computation using 
one-electron orbital wavefunctions for the core or target ion from an updated version of the 
atomic structure code SUPERSTRUCTURE, using a scaled Thomas- Fermi- Dirac potential 
(Eissner et al. 1974). 

The wavefunction expansion for O VI consists of 17 fine structure levels of configurations 
Is^, ls2s, ls2p, ls3s, ls3p, and ls3d of O VII. The levels, along with their relative energies, 
are given in Table 1. The set of correlation configurations in the atomic structure 
calculations and the Thomas- Fermi scaling paramenter (A„;) for each orbital are also given 
in Table 1. The second term in Eq. (3), with bound state correlation functions for O VI, 
includes all possible {N + l)-particle configurations with to maximum orbital occupancies: 
2s^ 2p2, 3s^ 3p^ 3d^ 4s and 4p. 

For O VII, the wavefunction expansion consists of 16 fine structure levels of O VIII of 
Is, 2s, 2p, 3s, 3p, 3d, 4s, 4p, Ad and 4/, as given in Table 1. The bound state correlation 
functions included all configurations from to maximum orbital occupancies: Is^, 2s^, 2p^, 
3s^, 3p^, 3d^, 4s^ and 4p^, 4d and 4/. The energies in Table 1 are observed values (from the 
NIST website: www.nist.gov). Although calculated energies are less than one percent of the 
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observed ones, the latter are used in the computations to obtain more accurate positions of 
resonances. Radial integrals for the partial wave expansion in Eq. 3 are specified for orbitals 
< £ < 9, with a R-matrix basis set 40 'continuum' functions (NRANG2) for O VI, and 30 
for O VII. 

Both the partial and the total photoionization cross setions are obtained for all bound 
levels. Coupled channel calculations for crpi include both the background and the resonance 
structures (due to the doubly excited autoionizing states) in the cross sections. Radiation 
damping of resonances up to n = 10 are included through use of the extended codes STGF 
and STGBF ( Nahar and Pradhan 1994, Zhang et al. 1999). The BPRM calculations are 
carried out for each total angular momentum symmetry Jtt, corresponding to a set of fine 
structure target levels Jj. Radiation damping of resonances within the close coupling BPRM 
calculations are described in Zhang et al. (1999, and references therein). The program 
PBRRAD is used to extend the total photoionization cross sections in the high energy 
region, beyond the highest target threshold in the close coupling wavefunction expansion 
of the ion, by a 'tail' using Kramers formula api{E) = aPI°{E°/EY where E° is the last 
tabulated energy above all target thresholds. 

Level-specific recombination cross sections crRc(^)) into bound levels i = n (SLJ) of 
the recombined (e + ion) system, are obtained from partial photoionization cross sections 
api{i,g) of the level i into the ground level g of the recombining ion. These detailed 
photo-recombination cross sections are calculated in the energy region from the threshold 
energy up to E^u = i/^ pa 10.0), where u is the effective quantum number relative to the 
target level of the recombining ion. The resonances up to u < i/g are delineated with 
a fine energy mesh. The electrons in this energy range generally recombine to a large 
number of final (e + ion) levels. Recombination cross sections are computed for all coupled 
symmetries and levels, and summed to obtain the total ctrc- Level specific recombination 
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rate coefficients are obtained using a new computer program, BPRRC (Nahar et al 2000). 
The level specific rates are obtained for energies going up to infinity. These rates include 
both non-resonant and resonant contributions up to energies Jul] Contributions from all 
autoionizing resonances up to < i^o ~ 10 are included. 

In the higher energy region, i/^ < < oo, below each target threshold where the 
resonances are narrow and dense and the background is negligible, we compute detailed 
and resonance averaged DR cross sections. The BPRM DR collision strengths are obtained 
using extensions of the i?- matrix asymptotic region code STGF (Nahar & Pradhan 1994, 
Zhang et al. 1999). It is necessary to use a very fine energy mesh in order to delineate the 
resonance structures. 

The program BPRRC sums up the level specific rates, which is added to the 
contributions from the resonant high-n DR, from resonances with < < oo, to obtain 
total recombination rates. As an additional check on the numerical calculations, the 
total recombination rate coefficients, aji, are also calculated from the total recombination 
collision strength, ^rc, obtained from all the photoionization cross sections, and the DR 
collision strengths. The agreement between the two numerical approaches is within a few 
percent. 

Finally, the small background (non-resonant) contribution from the high-n states 
(10 < n < oo) to total recombination is also included as the "top-up" part, computed 
in the hydrogenic approximation (Nahar 1996). This contribution is important at low 
temperatures, but negligible at high temperatures. At low temperature, recombination rate 
is dominated by the RR into the infinite number of high-n states, as electron energies are 
not usually high enough for resonant excitations and DR, and causes a rapid rise in q;r(T) 
towards low temperatures. 



4. RESULTS AND DISCUSSION 



Results are presented for photoionization and recombination of O VIII + e — > O VI, 
and O VIII + e — > O VII. Odd and even parity energy levels n (SLJ) of total angular 
momentum symmetries 1/2 < J < 17/2 are considered for O VI, and < J < 4 for O VII, 
all bound levels with n < 10. 

Photoionization cross sections api include both the total cross section, leaving the core 
(residual) ion in the ground and excited levels, as well as the partial cross sections, leaving 
the core ion in the ground level only. The total cross sections are needed in astrophysical 
applications, such as in ionization balance calculations while parital cross sections are 
needed for applications such as for recombination rate coefficients. (Both the total and the 
partial cross sections are available electronically.) 

Level-specific and total recombination rate coefficients using the BPRM unified 
treatment, (n SLJ, n < 10) and aRc(T) respectively, are presented for O VI and 
O VII. In the following subsections we describe salient features pertaining to aRc{T) in 
numerical form and in figures. Total aRc{T) are computed two different ways to enable 
numerical checks: (i) from the sum of the level-specific rate coefficients and the high-n DR 
contribution, and (ii) from total collision strengths calculated from photoionization cross 
sections directly, and the DR contribution. The differences between the two are typically 
within a few percent, thus providing a numerical and self-consistency check particularly on 
the resolution of resonances. We also present total recombination rate coefficients for the 
hydrogenic ion O VIII for completeness. 



4.1. O VI 



4.1.1. Photoionization 
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Figs. 1(a) and (b) show the ground state photoionization cross section for O VI 
(Is^ 2s ^5*1/2) • The top panel (a) presents the total photoionization cross section summed 
over the various target thresholds for ionization, and the bottom panel (b) presents the 
partial cross sections of the ground level into the ground Is^(^S') level of residual ion O VII. 
The resonances at high energies belong to Rydberg series converging n = 2, 3 levels. These 
are the well known KLL, KLn {n > 2) as, for example, discussed by Pradhan (2000) and 
Nahar et al. (2001). Since the first excited levels of n=2 thresholds of the core ion O VII lie 
at high energies, the cross sections decrease monotonically over a large energy range before 
the Rydberg series of resonances appears. The total and the partial cross sections are 
identical below the first excited level of the residual ion beyond which total apj increases 
due to added contributions from excited channels, as shown in the figures. A distinct 
difference between the total and partial cross sections Fig. 1 comes from the contribution 
of channels with excited n=2 thresholds. The K-shell ionization jump at the n = 2 target 
levels in total api is due to inner-shell photoionization: 

hu + VI{ls'^ 2s) — >e +0 VII{ls2s, ls2p) . 

In X-ray photoionization models inner-shell edges play an important role in overall 
ionization rates. 

Fig. 2 presents partial photoionization cross sections of the Rydberg series of levels 
Is^ns, 2 < n < 10. As noted in earher works, photoionization cross section of Rydberg 
levels exhibit wide resonances known as photoexcitation-of-core (PEC) resonances at 
energies associated with dipole transitions in the core ion. We see such PEC resonances in 
in cross sections of all bound levels of O VI, as in Fig. 2 at photon energies 41.788, 42.184, 
48.804 and 48.922 Ry due to core excitations to levels ls2p(3pf), ls2p(^Pf), ls3p(^Pf), and 
ls3p{^P°) of O VII. At these energies the core ion goes through an allowed transition, while 
the outer electron remains a 'spectator' in a doubly-excited resonance state, followed by 
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autoionization of the O VI resonance into the ground level of O VII. The effect gets more 
distinct for cross sections of higher excited levels. For low charge ions, PEC resonances are 
usualy wider. These resonances in O VI depict the behavior of cross sections that determine 
corresponding features in level-specific recombination rates. 

4.1.2. Recombination cross sections and rate coefficients 

The results consist of the total and 98 level-specific (n < 10) recombination rate 
coefficients q;r(T) and ai{T) respectively. Total recombination rates are given in Table 2. 
The main features are illustrated and compared with previously available data for RR and 
DR rates in Fig. 3. The sohd curve is the BPRM total unified aii{T) and shows typical 
features. The high values at very low temperature are due to the dominance of RR into 
an infinity of high-n levels. aR{T) decreases with increasing T until high temperatures 
where it rises due to the dominance of DR, followed by a monotonous decrease. The dotted 
curve in Fig. 3 is the total rate in LS couplig (Nahar 1999). The difference between the 
LS coupling aji{T) and the present BPRM rate is due to relativistic effects and radiation 
damping of low-n resonances. Both effects can be important for highly charged ions. The 
dashed curve is the RR rate by Pequignot et al (1991) obtained by smoothing of OP data 
over resonances, and hydrogcnic approximation. The agreement is good in the low region 
since non-resonant cross sections dominate the low energy region. The dot-dashed curve 
is the DR rate by Badnell et al. (1990) using the isolated resonance approximation. The 
sum of previous RR and DR rates agrees with the present unified rates to about 10-20%, as 
expected for highly charged ions. 

The 98 level-specific recombination rate coefficients ai{T) for O VI are for i = n{SLJ), 
n < 10 and i < 9, and associated Jtt levels. Fig. 4 presents ai{T) into eight n = 2 and 3 
levels with n (SLJ): 2s ^^o, 2p 3* ^-^o, 3p and 3d ^-03/2,5/2- These rates 
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are relatively smooth except for a small and diffuse DR 'bump'. The levels correspond to 
the prominent 1032 and 1038 A UV doublet hues respectively in O VI spectra arising due 
to strong dipole transitions ^-^3/2,1/2 — 'S'i/2- 

Total photorecombination cross sections a^c show features similar to photoionization 
cross sections api — smooth decay with energy, before the resonance complexes begin to 
emerge at high energies. However, {aRc are more complicated than apr since all levels 
contribute to the total auc, shown in Fig 5(a) for O VI. The resonance complexes are marked 
as KLL, KLM, KLN etc. For heavier high-Z elements than oxygen (Z > 10) these manifest 
themselves as di-electronic satellite lines observed in tokamaks, Electron-Beam-Ion- Traps 
(EBIT), ion storage rings and astrophysical sources. The KLL complexes have been well 
studied in previous works (e.g. Gabriel 1972, Bely-Dubau et al. 1982, Pradhan and Zhang 
1997, Zhang et al. 1999), Beiersdorfer et al. 1992, Oelgoetz and Pradhan 2001), for various 
ions. 

4.2. O VII 

4-2.1. Photoionization 

Total and partial photoionization cross sections apj for 116 levels n (SLJ) ( n < 
10) are computed (the partial cross sections are for ionization into the ground level 
ls(^5'i/2) of O VIII). Illustrative results are presented in Fig. 6. Fig. 6(a) shows the 
photoionization cross section of the ground level (Is^ ^5*0) of O VII. Similar to O VI, 
the Rydberg series of resonances, KL and Kn (n > 2), begins at fairly high energies 
owing to the high n = 2 excitation thresholds of O VIII. However, the ground level 
api of O VII does not show a significant K-shell jump at n = 2 threshold, as seen in 
O VI. Nonetheless the feature is prominent in the excited n — 2 level cross sections 
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as seen in Figs. 6(b)-(e). These excited levels are responsible for the formation of the 
important X-ray diagnostic lines, w, x, y, z, from the allowed and forbidden transitions 
ls2 {^So) i — Is2p(ipf),ls2p(3p°),ls2p(3pf),ls2s(35i), respectively. The K-shell 
ionization jump at the n = 2 target levels 

hu + O VII{ls2s, ls2p) — >e + VIII(2s , 2p) . 

can be seen clearly in the photoionization cross sections. 

4-2.2. Recombination cross sections and rate coefficients 

Total recombination rate coefficients for O VII over a wide temperature range are 
presented in Table 3. In Fig. 7 the solid curve is the present total recombination rate 
coefficients, including relativistic effects in BPRM approximation and radiation damping of 
low-n (n < 10) resonances. The sohd curve has the same general shape typical of unified 
recombination rates - the high rate at low temperature deceases until the recombination is 
dominated by DR forming a hump at high temperature. The dotted curve in the figure is 
the earlier total unified rate (Nahar 1999) where a part of the rates, the resonances 2lnl 
converging on to the n — 2 thresholds of O VIII, were also computed using the BPRM 
method including radiation damping (Zhang et al. 1999); the remainder was in LS coupling. 
Both results show very close agreement since the main n — 2 resonant contributions are 
the same. However, present rates also include n = 3, 4 resonant contributions Mnl, Aini. 
The dashed curve is the RR rate obtained by Verner and Ferland (1996), which agrees well 
with the present one in the region where DR is small. The dot-dashed curve is the DR 
rate calculated by Savin (1999) from experimentally measured recombination cross sections 
by Wolf et al. (1991). As is discernible from Fig. 7, the sum of the previous RR and 
the experimental DR rates (RR-I-DR) is significantly lower than the present unified rate 
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for O VII. At the temperature of peak DR contribution Log T = 6.6, the unified rate is 
about 20% higher the RR+DR value. This is most hkely because the experimental DR 
measurements cover a smaller energy range and do not include all contributions. 

Level-specific recombination rate coefficients are obtained for 116 levels n (SLJ) with 
< J < 4 and n < 10. These rates for fine structure levels were not obtained in our earlier 
work (Nahar 1999). Fig. 8 presents level specific rates for the n — 2 levels corresponding 
to the X-ray w, x, y, and z fines of O VII. The rates show a relatively smooth decay with 
temperature except for the high temperature DR bump for ls2p{^P^) and ls2s{^Si) levels. 

The total recombination cross sections anc, summed over recombination cross sections 
of individual fine structure levels up to n = 10 of O Vll and high-n contributions arc 
presented in Fig. 5(b). The resonance complexes LL, LM, LN, MM etc. of n = 2 and 3 
thresholds of the core are specified in the figure. The n — A resonances are too low for any 
significant contributions. These complexes are seen at high energies near the high lying 
O VIII thresholds. 

As mentioned above, as a numerical check on the calculations we verify that the sum of 
the level-specific rate coefficients and the DR contribution agrees within a few percent with 
the total recombination rate coefficient obtained from total coUision strengths obtained 
from cTflc in Fig. 5 for both O VI and O VII. 

4-2.3. Ionization fractions in coronal equilibrium 

Using the present total unified recombination rates, aji{T), for O VI and O VII we 
calculate ionization fractions of oxygen ions in coronal (coUisional) equilibrium, 

N{z - l,g)S{z - 1,5) = N{z,g)a{z,g), l<z< Zmax, (H) 
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where S is the rate coefficient for electron impact ionization and total element density 
Nt = Y.l='o N{z,g). The ionization fractions in -log[N{z) /Nt] are given in Table 3. 
However, present fractions show very little difference from earlier values (Nahar 1999). 
Most of the differencs are in the third significant digit and do not show up in Table 3. 

5. CONCLUSION 

Extensive results from relativistic calculations for total and level-specific photoionization 
and recombination cross sections and rates are presented for O VI and O VII. These are 
of general interest in UV and X-ray spectroscopy of laboratory and astrophysical sources, 
especially the formation due to recombination of important lines such as the AA 1032,1038 
fine structure doublet in O VI, and the 'triplet' X-ray features in O VII. 

A discussion of some of the important atomic effects such as resolution and radiation 
damping of resonances, interference between non-resonant (RR) and resonant (DR) 
recombination, comparison with experimental data and uncertainties, and general features 
of the unified (e -|- ion) recombination rates, has also been given in the first paper on the 
new BPRM calculations - paper IV in the present series (Nahar et al. 2000). 

Di-electronic satellite rates for the KLL, KLM, etc, complexes of several ions have 
earlier been shown to be in very good agreement with experiments and other theoretical 
calculations (Pradhan and Zhang 1997; Zhang et al. 1999), to about 10-20%; it is therefore 
expected that the present rates should be definitive, with similar uncertainties. 

The present level-specific data can be used to construct recombination-cascade matrices 
for O VI and O VII, to obtain effective recombination rates into specific fine structure levels 
n (SLJ) with n < 10 and £ < n — 1 (e.g. Pradhan 1985). The present data is more than 
sufficient for extrapolation to high-n,£ necessary to account for all cascade contributions. 
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The available data includes: (A) Photoionization cross sections, both total and partial, 
for bound fine structure levels of O VI and O VII up to the n = 10 levels. (B) Total 
unified recombination rates for O VI and O VII, and level-specific recombination rate 
coefficients for levels up to n = 10. Further calculations for other He-like and Li-like ions 
are in progress. All photoionization and recombination data are available electronically 
at: nahar@astronomy.ohio-state.edu. The total recombination rate coefficients are also 
available from the Ohio State Atomic Astrophysics website at: www.astronomy.ohio- 
state.edu/~pradhan. 

This work was supported partially by NSF and NASA. The computational work was 
carried out on the Cray SVl at the Ohio Supercomputer Center in Columbus, Ohio. 
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Table 1: Target levels in the eigenfunction expansions of O VII and O VIII. 





O VI 






O VIII 






level 


Et{Ry) 




level 


Et{Ry) 


1 


1 2/lo ^ 
Is^(^So) 


0.0 


1 


1 /2 O N 

isrSi/2) 


0.00 


2 


ls2s(^Si) 


41.232 


2 


2p(2p°/2) 


48.0308997 


3 


ls2p(:^P°(,) 


41.787 


3 


2s{2Si/2) 


48.0315698 


4 


Is2p(-'P°) 


41.788 


4 


2p('P3/2) 


48.0445917 


5 


ls2p(3p°) 


41.793 


5 


3p('P?/2) 


56.9304970 


6 


ls2s(lSo) 


41.812 


6 


3s(2Si/2) 


56.9306973 


7 


ls2p(lpf) 


42.184 


7 


3d(2Di/2) 


56.9345473 


8 


ls3s(3Si) 


48.651 


8 


3P('P3/2) 


56.9345541 


9 


ls3p(3Po) 


48.804 


9 


3d(2D5/2) 


56.9358976 


10 


Is3p(-^P°) 


48.804 


10 


4p('Pi/2) 


60.0447993 


11 


ls3p(^P^) 


48.804 


11 


4s(2Si/2) 


60.0448841 


12 


ls3s(iSo) 


48.811 


12 


4d(2D3/2) 


60.0465078 


13 


ls3d(3D3) 


48.884 


13 


4p{'P?/2) 


60.0465108 


14 


ls3d(3D2) 


48.884 


14 


4f('P?/2) 


60.0470765 


15 


ls3d(3Di) 


48.884 


15 


4d(2D5/2) 


60.0470775 


16 


ls3d(iDi) 


48.894 


16 


4f('P?/2) 


60.0473612 


17 


ls3p(ipj) 


48.922 









O VII: Correlations - 2s^ , 2p^ , 3s'^, 3p^, 3d^ , 2s2p, 2s3s, 

2s3p,2s3d, 2s4s, 2s4p, 2p3s, 2p3p, 2p3d, 2p4s, 2p4p, 

\„l - 0.991(ls),0.991(2s),0.776(2p),1.16883(3s),0.91077(3p), 

1.00746(3d),-1.59699(4s),-l-61237(4p) — see text) 

O VIII: A„; - 1.0, for Is to 4f 



- 21 - 



Table 2: Total recombination rate coefficients q;h(T) for O VI, O VII and O VIII. 





() \1 


() \ 11 


O \'iU 


(K) 


O \1 


() \11 


() Mil 


1.0 


1.42E- 


■09 


1.99E-09 


2.90E- 


■09 


5.1 


3.72E- 


■12 


6.16E- 


■12 


1.03E- 


■11 


1.1 


1.25E- 


■09 


1.76E-09 


2.56E- 


■09 


5.2 


3.13E- 


■12 


5.28E- 


■12 


8.78E- 


■12 


1.2 


I.IOE- 


■09 


1.55E-09 


2.26E- 


09 


5.3 


2.62E- 


■12 


4.51E- 


■12 


7.51E- 


■12 


1.3 


9.68E- 


■10 


1.36E-09 


2.00E- 


■09 


5.4 


2.19E- 


■12 


3.85E- 


■12 


6.42E- 


■12 


1.4 


8.50E- 


■10 


1.20E-09 


1.76E- 


■09 


5.5 


1.83E- 


■12 


3.28E- 


■12 


5.48E- 


■12 


1.5 


7.46E- 


■10 


1.05E-09 


1.55E- 


■09 


5.6 


1.52E- 


■12 


2.79E- 


■12 


4.67E- 


■12 


1.6 


6.53E- 


■10 


9.23E-10 


1.36E- 


■09 


5.7 


1.26E- 


■12 


2.37E- 


■12 


3.97E- 


■12 


1.7 


5.72E- 


■10 


8.10E-10 


1.20E- 


■09 


5.8 


1.05E- 


■12 


2.01E- 


■12 


3.37E- 


■12 


1.8 


5.01E- 


■10 


7.10E-10 


1.05E- 


■09 


5.9 


8.99E- 


■13 


1.72E- 


■12 


2.85E- 


■12 


1.9 


4.38E- 


■10 


6.21E-10 


9.26E- 


10 


6.0 


8.42E- 


13 


1.52E- 


■12 


2.40E- 


■12 


2.0 


3.83E- 


■10 


5.44E-10 


8.12E- 


■10 


6.1 


9.14E- 


■13 


1.44E- 


■12 


2.03E- 


■12 


2.1 


3.34E- 


■10 


4.75E-10 


7.12E- 


■10 


6.2 


l.llE- 


■12 


1.52E- 


■12 


1.71E- 


■12 


2.2 


2.91E- 


■10 


4.15E-10 


6.23E- 


■10 


6.3 


1.39E- 


■12 


1.73E- 


■12 


1.43E- 


■12 


2.3 


2.54E- 


■10 


3.62E-10 


5.46E- 


■10 


6.4 


1.65E- 


■12 


1.99E- 


■12 


1.20E- 


■12 


2.4 


2.21E- 


■10 


3.16E-10 


4.78E- 


■10 


6.5 


1.82E- 


■12 


2.21E- 


■12 


9.96E- 


13 


2.5 


1.93E- 


■10 


2.76E-10 


4.18E- 


■10 


6.6 


1.86E- 


■12 


2.31E- 


■12 


8.27E- 


13 


2.6 


1.67E- 


■10 


2.40E-10 


3.65E- 


■10 


6.7 


1.78E- 


■12 


2.26E- 


■12 


6.84E- 


13 


2.7 


1.46E- 


■10 


2.09E-10 


3.19E- 


■10 


6.8 


1.60E- 


■12 


2.10E- 


■12 


5.63E- 


■13 


2.8 


1.27E- 


■10 


1.82E-10 


2.78E- 


■10 


6.9 


1.38E- 


■12 


1.85E- 


■12 


4.61E- 


■13 


2.9 


I.IOE- 


■10 


1.58E-10 


2.43E- 


■10 


7.0 


1.14E- 


■12 


1.56E- 


■12 


3.75E- 


■13 


3.0 


9.52E- 


■11 


1.37E-10 


2.12E- 


10 


7.1 


9.19E- 


13 


1.28E- 


■12 


3.05E- 


13 


3.1 


8.25E- 


■11 


1.19E-10 


1.85E- 


■10 


7.2 


7.20E- 


■13 


1.02E- 


■12 


2.45E- 


13 


3.2 


7.15E- 


■11 


1.03E-10 


1.61E- 


■10 


7.3 


5.54E- 


■13 


7.96E- 


13 


1.96E- 


13 


3.3 


6.19E- 


■11 


8.97E-11 


1.40E- 


■10 


7.4 


4.19E- 


13 


6.10E- 


13 


1.57E- 


■13 


3.4 


5.35E- 


■11 


7.77E-11 


1.22E- 


■10 


7.5 


3.13E- 


13 


4.61E- 


13 


1.24E- 


■13 


3.5 


4.62E- 


■11 


6.73E-11 


1.06E- 


■10 


7.6 


2.32E- 


13 


3.45E- 


13 


9.76E- 


■14 


3.6 


3.99E- 


■11 


5.82E-11 


9.22E- 


■11 


7.7 


1.70E- 


13 


2.56E- 


13 


7.63E- 


■14 


3.7 


3.43E- 


■11 


5.04E-11 


8.01E- 


■11 


7.8 


1.24E- 


■13 


1.89E- 


■13 


5.93E- 


■14 


3.8 


2.96E- 


11 


4.35E-11 


6.95E- 


11 


7.9 


9.00E- 


■14 


1.38E- 


13 


4.59E- 


■14 


3.9 


2.54E- 


■11 


3.76E-11 


6.03E- 


■11 


8.0 


6.50E- 


■14 


I.OIE- 


■13 


3.51E- 


■14 


4.0 


2.18E- 


■11 


3.25E-11 


5.22E- 


■11 


8.1 


4.69E- 


■14 


7.35E- 


■14 


2.70E- 


■14 


4.1 


1.87E- 


■11 


2.80E-11 


4.53E- 


■11 


8.2 


3.37E- 


■14 


5.33E- 


■14 


2.05E- 


■14 


4.2 


1.61E- 


■11 


2.41E-11 


3.92E- 


■11 


8.3 


2.42E- 


■14 


3.86E- 


■14 


1.55E- 


■14 


4.3 


1.37E- 


■11 


2.08E-11 


3.39E- 


■11 


8.4 


1.73E- 


■14 


2.79E- 


■14 


1.16E- 


■14 


4.4 


1.17E- 


■11 


1.79E-11 


2.93E- 


■11 


8.5 


1.24E- 


■14 


2.02E- 


■14 


8.72E- 


■15 


4.5 


l.OOE- 


■11 


1.54E-11 


2.53E- 


■11 


8.6 


8.87E- 


■15 


1.45E- 


■14 


6.51E- 


■15 


4.6 


8.54E- 


■12 


1.32E-11 


2.18E- 


■11 


8.7 


6.33E- 


■15 


1.05E- 


■14 


4.83E- 


■15 


4.7 


7.26E- 


■12 


1.14E-11 


1.88E- 


■11 


8.8 


4.52E- 


■15 


7.55E- 


■15 


3.58E- 


■15 


4.8 


6.16E- 


■12 


9.77E-12 


1.62E- 


■11 


8.9 


3.23E- 


■15 


5.44E- 


■15 


2.64E- 


■15 


4.9 


5.22E- 


■12 


8.38E-12 


1.39E- 


11 


9.0 


2.30E- 


■15 


3.91E- 


■15 


1.93E- 


■15 


5.0 


4.41E- 


■12 


7.20E-12 


1.19E- 


■11 
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Table 3: Ionization fractions, —logio^^, of oxygen ions in coronal equilibrium. 





O I 


O II 


O III 


O IV 


logioN(z)/NT 
O V O VI 


O VII 


O VIII 


O IX 


4.0 


0.001 


2.5 


***** 


***** 


***** 


***** 


***** 


***** 


***** 


4.1 


0.0433 


1.02 


***** 


***** 


***** 


***** 


***** 


***** 


***** 


4.2 


0.432 


0.2 


7.91 


***** 


***** 


***** 


***** 


***** 


***** 


4.3 


1.19 


0.029 


5.37 


***** 


***** 


***** 


***** 


***** 


***** 


4.4 


1.77 


0.008 


3.49 


***** 


***** 


***** 


***** 


***** 


***** 


4.5 


2.12 


0.007 


2.10 


8.15 


***** 


***** 


***** 


***** 


***** 


4.6 


2.38 


0.037 


1.11 


5.38 


***** 


***** 


***** 


***** 


***** 


4.7 


2.72 


0.176 


0.481 


3.39 


9.73 


***** 


***** 


***** 


***** 


4.8 


3.21 


0.483 


0.179 


2.05 


6.89 


***** 


***** 


***** 


***** 


4.9 


3.82 


0.901 


0.096 


1.14 


4.76 


***** 


***** 


***** 


***** 


5.0 


4.52 


1.40 


0.174 


0.539 


3.17 


7.54 


***** 


***** 


***** 


5.1 


5.33 


2.01 


0.417 


0.23 


2.02 


5.13 


8.58 


***** 


***** 


5.2 


6.22 


2.71 


0.778 


0.114 


1.21 


3.28 


5.48 


***** 


***** 


5.3 


7.16 


3.45 


1.22 


0.155 


0.646 


1.87 


3.04 


***** 


***** 


5.4 


8.23 


4.32 


1.82 


0.41 


0.387 


0.895 


1.23 


***** 


***** 


5.5 


9.78 


5.67 


2.92 


1.20 


0.724 


0.634 


0.289 


***** 


***** 


5.6 


3|e 3|e :4: ^ 


7.43 


4.43 


2.43 


1.57 


0.975 


0.064 


8.40 


***** 


5.7 


***** 


9.11 


5.89 


3.64 


2.43 


1.40 


0.019 


6.32 


***** 


5.8 


***** 


***** 


7.19 


4.71 


3.18 


1.78 


0.007 


4.66 


***** 


5.9 


***** 


***** 


8.33 


5.63 


3.82 


2.10 


0.004 


3.31 


7.86 


6.0 


***** 


***** 


9.31 


6.41 


4.34 


2.34 


0.005 


2.23 


5.52 


6.1 


***** 


***** 


***** 


7.05 


4.75 


2.49 


0.020 


1.39 


3.66 


6.2 


***** 


***** 


***** 


7.61 


5.09 


2.60 


0.080 


0.795 


2.24 


6.3 


***** 


***** 


***** 


8.19 


5.47 


2.76 


0.231 


0.455 


1.22 


6.4 


***** 


***** 


***** 


8.89 


5.98 


3.08 


0.527 


0.368 


0.56 


6.5 


***** 


***** 


***** 


9.74 


6.65 


3.56 


0.984 


0.51 


0.231 


6.6 


***** 


***** 


***** 


***** 


7.41 


4.15 


1.52 


0.77 


0.096 


6.7 










8.16 


4.75 


2.05 


1.06 


0.044 


6.8 










8.90 


5.33 


2.55 


1.33 


0.022 


6.9 










9.59 


5.88 


3.00 


1.57 


0.012 


7.0 












6.41 


3.42 


1.79 


0.007 


7.1 












6.91 


3.81 


1.98 


0.005 


7.2 












7.39 


4.17 


2.16 


0.003 


7.3 












7.85 


4.50 


2.32 


0.002 


7.4 












8.28 


4.82 


2.47 


0.001 


7.5 












8.70 


5.12 


2.61 


0.001 
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lO 15 20 25 30 35 40 45 50 55 

Photon Energy (Ry) 



Fig. 1. — Photoionization cross sections of the ground level Is^ 2s of O VI: (a) Total 

cross section; the large jump around 52 Ry corresponds to the n=2 K-shell ionization edge, 
(b) Partial cross section into the ground level Is^ (^-S'o) of O VII; note that the jump is no 
longer present and the cross section is continuous across the n = 2 thresholds of O VII. 
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Fig. 2. — Partial photoionization cross sections of the Rydberg series of levels ls^ns(^Si/2) of 
O VI into the ground state state 1s'^{^Sq) of O VII. photoexcitation-of-core (PEC) resonances 
are seen at about 42 and 49 Ry. 
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Fig. 3. — Total unified recombination rate coefficients for O VI: BPRM with fine structure 
(sofid curve); Total in LS coupling (Nahar 1999, dotted): RR rates by Pequignot et al (1991, 
dash), DR rates by Badnell et al. (1990, dot-dash). 
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2 4 6 8 

logioT(K) 



Fig. 4. — Level-specific recombination rate coefficients for O VI recombining to ground and 
excited n—2, 3 levels. The strong dipole transitions ^^3° 2 1/2 — ^1/2 responsible for 
the prominent 1032 and 1038 A UV doublet lines respectively in O VI spectra. 
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Fig. 5. — Total unified (e + ion) photo-recombination cross sections, grc^ of (a) O VI 
and (b) O VII. Note that the grc exhibit considerably more resonance structures than the 
corresponding ground level cxp/ in Figs. 1 and 6, since the former are summed over the 
ground and many excited recombined levels. 



_l 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

O VII + hv -> O vm + e 


1 1 1 1 i_ 


= (a) Is^('So): Total 

H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 










1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


1 n=2 


1 1 1 1 r 



lO 

1 
.1 

.Ol 

.OOl 
.OOOl 

1 
.1 

.Ol 
.OOl 
.OOOl 

.1 

.Ol 
.OOl 
.OOffiffi 

1 
.1 

.Ol 
.OOl 
.OOOl 



55 60 6.5 7P 75 8p 85 90 



i ^r-^ i I i !V? i i i 1V-7 i i llV i 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

15 20 25 30 35 40 45 50 55 60 65 70 75 

Photon Energy (Ry) 



Fig. 6. — Level-specific photoionization cross sections of (a) the ground level Is^ (^•S'o), and 
(b) - (e) excited ls2s and ls2p levels of O VII responsible for the prominent X-ray lines: 
resonance (w), intercombination (x,y), and forbidden(z). 
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Fig. 7. — Total unified recombination rate coefficients for O VII: present BPRM (solid); 
total (dotted, Nahar 1999), RR rates by Verner and Ferland (1996) (dahsed), DR rates by 
Savin (1999) from measured cross sections by Wolf et al. (1991) (dot-dashed). 
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Fig. 8. — Level-specific recombination rate coefficients for O VII into the excited n — 2 levels 
responsible for the prominent X-ray lines w, x, y, and z. 



